2D symmetry protected topological orders and their protected gapless edge excitations 



Xie Chen, 1 Zheng-Xin Liu, 1,2 and Xiao-Gang Wen 1 ' 2 

1 Department of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA 
2 Institute for Advanced Study, Tsinghua University, Beijing, 100084, P- R- China 

Topological insulators in free fermion systems have been well characterized and classified. How- 
ever, it is not clear in strongly interacting boson or fermion systems what symmetry protected 
topological orders exist. In this paper, we present a model in a 2D interacting spin system with 
nontrivial on-site Z2 symmetry protected topological order. The order is nontrivial because we can 
prove that the ID system on the boundary must be gapless if the symmetry is not broken, which 
generalizes the gaplessness of Wess-Zumino-Witten model for Lie symmetry groups to any discrete 
symmetry groups. The construction of this model is related to a nontrivial 3-cocycle of the Z2 group 
and can be generalized to any symmetry group. It potentially leads to a complete classification of 
symmetry protected topological orders in interacting boson and fermion systems of any dimension. 
Specifically, this exactly solvable model has a unique gapped ground state on any closed manifold 
and gapless excitations on the boundary if Z2 symmetry is not broken. We prove the latter by devel- 
oping the tool of matrix product unitary operator to study the nonlocal symmetry transformation 
on the boundary and revealing the nontrivial 3-cocycle structure of this transformation. Similar 
ideas are used to construct a 2D fermionic model with on-site Z2 symmetry protected topological 
order. 

PACS numbers: 71.27.+a, 02.40.Re 



I. INTRODUCTION 

Topological phases of matter are gapped quantum sys- 
tems containing nontrivial orders which are not due to 
spontaneous symmetry breaking in the ground states. 
While topologically ordered systems all have exponen- 
tially decaying correlation and appear quite simple from 
a classical point of view, various exotic quantum features 
have been discovered which reveal the surprisingly rich 
structure of topologically ordered systems. For example, 
some systems have ground state degeneracy which de- 
pends on the topology of the closed manifold the system 
is on; 1 ' 2 some have protected gapless edge excitations if 
the system has a boundary; 3,4 some have nontrivial en- 
tanglement structure in the ground state; 5,6 and some 
have bulk excitations with nontrivial statistics. How to 
obtain a clear picture of topological phases among such 
a variety of phenomena? First we find that topological 
phases can be divided into two general classes according 
to its level of stability under perturbations. 

The first class has 'intrinsic' topological order. 1 Sys- 
tems in this class must go through a phase transition to 
a trivial phase no matter what kind of local perturba- 
tion is added. Or using the local unitary equivalence 
between ground states we find that this class of sys- 
tems have ground states which cannot be mapped to a 
product state under ANY local unitary transformation 
as defined in Ref. 10. We say that this kind of states 
have long range entanglement. Example systems in this 
class include quantum Hall(integer or fractional), 11,12 
p + ip superconductor, 13,1 string-net models, 9 Z 2 spin 
liquid, 15-17 and chiral spin liquid. 18,19 It has been discov- 
ered that systems with 'intrinsic' topological order usu- 
ally have topology dependent ground state degeneracy, 1 ' 2 
nontrivial topological entanglement entropy 5 ' 6 and frac- 



tional statistics of bulk excitation. 7-9 In the following dis- 
cussion we will use the term 'topological order' to specif- 
ically refer to this class of systems. 

The second class has 'symmetry protected' topological 
order. This kind of system has certain symmetry and 
its non-degenerate ground state does not break any of 
the symmetries. If arbitrary perturbations are allowed, 
systems in this class all belong to the same phase as a 
trivial state. Its ground state can be mapped to a product 
state with local unitary transformations and hence are 
short range entanglcd(SRE). However, if only symmetric 
perturbations are allowed, systems in this class are in 
different phases from the trivial phase. Therefore, we 
say that the topological order in this class is symmetry 
protected. We will call these phases 'symmetry protected 
topological' (SPT) phases. Example systems in this class 
include Haldane phase in one dimensional spin chain 20 
and topological insulators. 21-26 Systems with SPT order 
have non-degenerate ground states on closed manifold 
and usually have nontrivial edge degrees of freedom if 
the system has a boundary. 2 1-27 

Many efforts have been made to obtain a more com- 
plete understanding of topological and symmetry pro- 
tected topological orders. In particular, topological 
and SPT orders have been completely classified in one- 
dimensional spin systems. 28,29 It was found that one- 
dimensional spin systems cannot have nontrivial topo- 
logical order but different SPT orders exist for sys- 
tems with certain symmetry. Similarly, a classifica- 
tion of fermion systems (interacting) in one dimension 
is also possible. 30-32 The picture changes dramatically in 
higher dimensions. First of all, nontrivial topological or- 
der does exist in two or higher dimensions. A lot has 
been learned about possible topological orders 9 ' 10,33-35 
although a complete understanding is still missing. In 
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this paper, we are going to focus only on the SPT phases. 
Most SPT phases in two and higher dimensions have 
been identified in free fcrmion systems due to the sim- 
plicity and versatility of the formalism. A classification 
of possible SPT phases in non-interacting fermion sys- 
tems has been obtained. 36-38 The major open question 
about SPT phases is in general which of these phases re- 
main and what new SPT phases are possible when the 
system is strongly interacting. In boson systems, even 
less is known as non-interacting bosons are necessarily 
topologically trivial. 39 

In this paper, we present a generic picture for under- 
standing SPT phases in interacting systems through the 
explicit construction of a simple example. Instead of 
starting from free fermions, we take a different approach 
and generalize our understanding of one dimensional in- 
teracting SPT phases to construct a two dimensional spin 
model with on-site Zi symmetry protected topological 
order. We call this model the CZX model for reason 
that will become clear later. On a closed surface the 
CZX model looks simple. Its Hamiltonian is composed 
of commuting projectors. Its symmetric gapped ground 
state is a product of local loops and hence short range 
entangled. However, the model becomes highly nontriv- 
ial if it has a boundary. The boundary must have gapless 
excitation as long as symmetry is not broken, a signature 
of nontrivial SPT order. We prove this fact by relating 
effective symmetry transformation on the boundary with 
a nontrivial 3-cocycle of the Z-i group. 

The construction of the CZX model signifies the close 
relation between SPT phases and nontrivial cocycles of 
the symmetry group. This idea is not limited to two 
dimensional systems. In another paper, 40 we generalize 
the formalism and construct nontrivial SPT phases in 
any d dimension with on-site unitary and anti-unitary 
symmetries G based on (d+ l)-cocycles of G. We expect 
that this construction gives a complete classification of 
d-dimensional SPT phases. 

The effective theory on the boundary can be seen 
as a generalization of the Wess-Zumino-Witten (WZW) 
model. 41 ' 42 The WZW model describes conformally in- 
variant ID systems with an internal symmetry of a com- 
pact Lie group. The WZW model obtained by adding 
a topological term (the WZW term) to the usual dy- 
namical term in the Lagrangian of the nonlinear sigma 
model, is exactly solvable in semiclassical limit. It ex- 
plains the physics of ID gapless systems with a global 
Lie group symmetry. However, the construction of the 
model depends crucially on the fact that the symmetry 
group is continuous and does not apply to, for exam- 
ple, the Z2 group. Our proof of the gapless-ness of the 
ID effective theory on the boundary of the CZX model 
hence generalizes the understanding of the WZW model 
to discrete groups. Our method based on the nontrivial 
3-cocycles applies to both continuous and discrete sym- 
metry groups, although it does not give the conformal 
field theory of the system directly. Also our proof is non- 
perturbative, not relying on semiclassical approximation. 



The connection between the CZX model and the WZW 
model is not particularly clear in the formulation of this 
paper, as the WZW model is usually given in the La- 
grangian form. In another paper 40 , we reformulate our 
models (including the CZX model and those for all other 
symmetries and in all dimensions) in the Lagrangian lan- 
guage where the connection with the WZW model would 
become obvious. 

The paper is organized as follows: in section II, we re- 
view our understanding of the entanglement structure of 
SPT phases in one dimension. In generalizing such entan- 
glement structure to higher dimension, we first present a 
naive attempt which fails to produce interesting phases. 
Identifying the missing element, we construct the CZX 
model in section III. We give explicitly the symmetry of 
the system, its Hamiltonian and its ground state. In or- 
der to show the nontrivial-ness of this model, we study 
its effective boundary theory in section IV. We identify 
the effective degrees of freedom, effective Zi symmetry 
and show that in simple cases the boundary cannot be in 
a gapped symmetric phase. In order to prove this conclu- 
sion in general, we use the tool of matrix product unitary 
operators(MPUO). Introduction to the matrix product 
unitary operators formalism is given in appendix D in- 
cluding its definition and some simple properties. In sec- 
tion V, we show how to represent the effective symmetry 
on the boundary of the CZX model using MPUO. We find 
that the transformation rule between the MPUO's is re- 
lated to a nontrivial class of 3-cocycles in the third coho- 
mology group % 3 (Z 2 , C/(l)) of Z 2 . 43 Using this relation, 
we prove that the boundary cannot have a gapped sym- 
metric ground state. This result applies in general to any 
MPUO related to a nontrivial 3-cocycle in n 3 (G, U(l)). 
Hence we conclude that the CZX model is in a nontriv- 
ial SPT phase protected by on-site Z 2 symmetry. Using 
similar ideas, we construct in section VII a fermion sys- 
tem with on-site Z2 symmetry whose boundary is also 
nontrivial. 



II. FROM ID SPT PHASES TO 2D 

In this section we first review our understanding of 
the entanglement pattern at the fixed point of ID SPT 
phases which we then try to generalize to higher dimen- 
sions. However, we are going to show that a straight 
forward generalization fails to give nontrivial SPT order. 
We identify the missing elements and prepare for the con- 
struction of nontrivial model in the next section. 

Each ID SPT phase in systems with on-site symme- 
try G can be well understood from the entanglement 
pattern of its ground state at fixed point, as shown in 
Fig.l. At fixed point, each site contains two spins. On 
each site, symmetry is represented linearly. But on each 
spin, symmetry only needs to be represented projectively. 
(A simple example of projective representation is given 
by SO(3) symmetry on a spin 1/2. For an introduction 
to projective representations and the second cohomol- 
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ogy group "H 2 (G, i7(l)) see appendix A. More generally, 
group cohomology is introduced in appendix B.) If sym- 
metry on the left spin belongs to the projective represen- 
tation of class u> in H 2 (G, £7(l))(for example spin 1/2 un- 
der 5*0(3)), then on the right spin it belongs to (again 
spin 1/2 under 50(3)) so that together they form a linear 
representation. The ground state of the system is a prod- 
uct of dimers between spins on neighboring sites. Each 
dimer is an entangled state of two spins which forms a 
one dimensional representation of G. The ground state is 
hence a total singlet under the symmetry. The nontrivial 
feature of the system shows up when we cut the chain 
into a finite segment. There are free degrees of freedom 
at the ends of the segment, each forming a projective 
representation of G. Two ID systems belong to the same 
SPT phase if their end degrees of freedom belong to the 
same class of projective representation uj g H 2 (G, £7(1)). 

CjO -CjJ UJ -U) U) -U) U) -U) 

C S^ dX) dj> (S > 

FIG. 1. Fixed point ground state of ID SPT phase with on- 
site symmetry of group G. Each site contains two spins, which 
form projective representation of class ui and —ui respectively. 
Connected spins form a dimer which forms a one-dimensional 
representation of G. On a finite segment of the ID chain, the 
boundary spins form projective representations of G. 

This simple picture can be generalized to two or higher 
dimension to give a 'bond' state. Consider the 2D state 
in Fig. 2. 




FIG. 2. A 2D 'bond' state which is short range entangled 
and is symmetric under on-site symmetry of group G. Each 
site contains four spins, each forming a projective represen- 
tation of G. Two spins connected by a bond form projective 
representations of class to and —u) respectively. The 'bond' 
represents an entangled state of the two spins which forms 
an one-dimensional representation of G. On a lattice with 
boundary, the boundary degrees of freedom are spins with 
projective representation uj(—uj.) 

Every site contains four spins. Each spin forms a pro- 
jective representation of on-site symmetry G, but the four 
spins on each site together form a linear representation of 
G. Two spins on neighboring sites which are connected 
by a bond forms projective representation uj and — to re- 
spectively and the bond represents an entangled state 



between the two spins which forms a one dimensional 
representation of G. Similar to the ID case, the total 
state is invariant under on-site symmetry G. The state 
is short range entangled and can be the gapped ground 
state of a simple Hamiltonian(sum of projections onto 
the entangled pairs). If the system is defined on a disk 
with boundary, there will be free degrees of freedom at 
each site on the boundary which form projective repre- 
sentations of G. 

It might seem that states with different projective rep- 
resentations at each site on the boundary correspond to 
different SPT phases, just like in the ID case. How- 
ever, this is not totally true. If translation symmetry is 
required, each boundary spin is well defined and the pro- 
jective representation they form do label different phases. 
On the other hand, in the absence of translation sym- 
metry, boundary spins can be combined and their pro- 
jective representations can add together. As projective 
representations form an additive group (the second coho- 
mology group "H 2 (G, £7(1)) of G), combining boundary 
spins would change the projective representations from 
one class to another and in particular, to the trivial class. 
Therefore, without translation symmetry, all 2D states 
with a bond form as shown in Fig. 2 belong to the same 
phase. 

On the other hand, SPT phases are known to exist 
in two and higher dimensions without the protection of 
translation symmetry, for example in topological insula- 
tors. The simple bond picture above therefore cannot 
account for their SPT order. In order to have nontriv- 
ial SPT order, we need to generalize the bond state in 
two ways: (1) the local entanglement structure is not 
bonds between two spins, but rather plaquettes among 
four spins on sites around a square. This alone is not 
enough to construct new SPT order. We also need (2) 
symmetry transformation on each site does not factorize 
into separate operations on each of the four spins. That 
is, the total linear symmetry operation on each site is 
not a tensor product of four projective representations 
as otherwise the state can be reduced to a bond state. 

Following this line of thought, we construct the CZX 
model in section III. The CZX model has an on-site Z-i 
symmetry that does not factorize into projective repre- 
sentations and the symmetry protected topological order 
of the state is robust against disorder. The boundary ef- 
fective degrees of freedom in CZX model has an effective 
Zi symmetry which cannot be written in an on-site form. 
Moreover, the boundary cannot be in a gapped symmet- 
ric state under the effective symmetry. In other words, 
the boundary must either break the Z% symmetry or have 
gapless excitations. This is different from the bond state 
discussed above(Fig.2). In the bond state, the boundary 
degrees of freedom are the boundary spins with projective 
representations. The effective symmetry is still on-site. 
Several boundary spins can form a singlet if their pro- 
jective representations add up to a linear representation. 
Therefore, in the bond state, the boundary can be in a 
gapped symmetric state under on-site symmetry simply 
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by breaking translation symmetry. However, in the CZX 
model, this is not possible. 

III. CZX MODEL 

In this section, we construct the CZX model explicitly 
which turns out to have nontrivial SPT order protected 
only by on-site Z 2 symmetry. 




FIG. 3. CZX model (a) each site (circle) contains four 
spins (dots) and the spins in the same plaquette (square) 
are entangled. (b) on-site Z 2 symmetry is generated by 
Uczx = X 1 X 2 X z X i CZ 12 CZ 2 3CZ zi CZ il (c) a local term 
in the Hamiltonian, which is a tensor product of one X4 term 
and four P 2 terms as defined in the main text. 

Consider a square lattice with four two-level spins per 
site, as shown in Fig. 3(a) where sites are represented by 
circles and spins are represented by dots. We denote the 
two levels as |0) and |1). The system has an on-site Z 2 
symmetry as given in Fig. 3(b). It is generated by 

Uczx = U X U CZ (1) 

where 

U x = X x ® X 2 ® X 3 ® X 4 (2) 

Xi is Pauli X operator on the zth spin and 

Ucz = CZ12CZ23CZ34CZ41 (3) 

where CZ is the controlled-.^ operator on two spins de- 
fined as 

CZ = |00)(00| + |01)(01| + |10)(10| - |11)(11| (4) 

As defined, CZ does nothing if at least one of the spins is 
in state |0) and it adds a minus sign if both spins are in 
state |1). Different CZ operators overlap with each other. 
But because they commute, Ucz is well defined. Note 
that Ucz cannot be decomposed into separate operations 
on the four spins and the same is true for Uczx- Ux and 
Ucz both square to / and they commute with each other. 
Therefore, Uczx generates a Z2 group. 



The Hamiltonian of the system is defined as a sum of 
local terms around each plaquette. Plaquettes are repre- 
sented by squares in Fig. 3. H = ^ H p ., where the term 
around the ith plaquette H Pi acts not only on the four 
spins in the plaquette but also on the eight spins in the 
four neighboring half plaquettes as shown in Fig. 3(c) 

H Vi = -Xi ® P 2 " ® F 2 d «> P l 2 ® P r 2 (5) 

where X4 acts on the four spins in the middle plaquette 

as 

X4, = |0000)(1111| + |1111)(0000| (6) 

and P2 acts on the two spins in every neighboring half 
plaquette as 

P 2 = |00)(00| + |11)<11| (7) 

P 2 , P 2 , Pi, P-2 acts on the up, down, left and right 
neighboring half plaquettes respectively. For the remain- 
ing four spins at the corner, H Pi acts as identity on them. 
The P2 factors ensure that each term in the Hamiltonian 
satisfies the on-site Z 2 symmetry defined before. 

All the local terms in the Hamiltonian commute with 
each other, therefore it is easy to solve for the ground 
state. If the system is defined on a closed surface, it has 
a unique ground state which is gapped. In the ground 
state, every four spins around a plaquette are entangled 
in the state 

\i> Pi ) = 10000) + |1111) (8) 

and the total wavefunction is a product of all plaquette 
wavefunction. If we allow any local unitary transforma- 
tion, it is easy to see that the ground state can be disen- 
tangled into a product state, just by disentangling each 
plaquette separately into individual spin states. There- 
fore, the ground state is short range entangled. However, 
no matter what local unitary transformations we apply 
to disentangle the plaquettes, they necessarily violate the 
on-site symmetry and in fact, the plaquettes cannot be 
disentangled if the Z 2 symmetry is preserved, due to the 
nontrivial SPT order of this model which we will show in 
the next sections. 

It can be checked that this ground state is indeed in- 
variant under the on-site Z 2 symmetry. Obviously this 
state is invariant under Ux applied to every site. It is 
also invariant under Ucz applied to every site. To see 
this note that between every two neighboring plaquettes, 
CZ is applied twice, at the two ends of the link along 
which they meet. Because the spins within each plaque- 
tte are perfectly correlated (they are all |0) or all |1)), 
the effect of the two CZ's cancel each other, leaving the 
total state invariant. 

Therefore, we have introduced a 2D model with on- 
site Z 2 symmetry whose ground state does not break the 
symmetry and is short-range entangled. In particular, 
this on-site symmetry is inseparable as discussed in the 
introduction and therefore cannot be characterized by 



5 



projective representation as in the bond state. We can 
add small perturbation to the system which satisfies the 
symmetry and the system is going to remain gapped and 
the ground state short range entangled and symmetric. 
It seems that the system is quite trivial and boring. How- 
ever, we are going to show that surprising things happen 
if the system has a boundary and because of these special 
features the system cannot be smoothly connected to a 
trivial phase even if translation symmetry is not required. 



IV. CZX MODEL BOUNDARY 

The non-trivialness of this model shows up on the 
boundary. Suppose that we take a simply connected disk 
from the lattice, as shown in Fig. 4(a). 




FIG. 4. (a)CZX model on a disk with boundary (b) boundary 
effective degrees of freedom form a ID chain which cannot 
have a SRE symmetric state (c) two boundaries together can 
have a SRE symmetric state which is a product of entangled 
pairs between effective spins connected by a dashed line. 

The reduced density matrix of spins in this region is 
invariant under on-site symmetry in this region. The 
reduced density matrix is a tensor product of individual 
terms on each full plaqucttc, half plaquette and corner of 
plaquette respectively. On a full plaquette 



P4 = (|oooo> + |mi))((oooo| + (iiii|) 



On a half plaquette 



p 2 = |00)(00| + |11)(11| 



On a corner of a plaquette 



Pi = 



umi 



(9) 



(10) 



(11) 



The state of spins on the plaqucttcs totally inside this 
region is completely fixed. But on the boundary there 
are free degrees of freedom. However, unlike in the bond 
state, only part of the total Hilbert space of the spins 
on the boundary is free. In particular, two spins in a 
half plaquette on the boundary are constrained to the 
two-dimensional subspace |00) (00 1 + |11)(11| and form 
an effective spin degree of freedom if we map 1 00) to |0) 
and |11) to |1). 

In Fig. 4(b), we show the effective degrees of freedom 
on the boundary as diamonds on a line. Projecting the 



total symmetry operation on the disk to the space sup- 
porting reduced density matrix, we find that the effec- 
tive symmetry operation on the boundary effective spins 
is Uczx = Jlili X i riili w ith Pauli X on each 

effect spin and CZ operation between neighboring effec- 
tive spins. The boundary is periodic and CZn,n+i acts 
on effective spin N and 1. This operator generates a Z 2 
symmetry group. 

This is a very special symmetry on a ID system. First 
it is not an on-site symmetry. In fact, no matter how 
we locally group sites and take projections, the sym- 
metry operations are not going to break down into an 
on-site form. Moreover, no matter what interactions we 
add to the boundary, as long as it preserves the sym- 
metry, the boundary cannot have a gapped symmetric 
ground state. We can start by considering some simple 
cases. The simplest interaction term preserving this sym- 
metry is ZiZ i+ i. This is an Ising interaction term and 
its ground state breaks the Z2 symmetry. In the trans- 
verse Ising model, the system goes to a symmetric phase 
if magnetic field in the x direction is increased. However, 
Xi breaks the Z 2 symmetry Uczx on the boundary and 
therefore cannot be added to the Hamiltonian. In fact, 
we are going to prove that the boundary cannot have 
SRE symmetric ground state (actually a more general- 
ized version of it) in the next section. This is one special 
property that differs the CZX model from the bond state 
in Fig. 2. In the bond state, the symmetry operations on 
the boundary are just projective representations on each 
site. Without translational invariance, there can always 
be a SRE symmetric state with this symmetry. 

The special property on the boundary only shows up 
when there is an isolated single boundary. If we put two 
such boundaries together and allow interactions between 
them, everything is back to normal. As shown in Fig. 4(c), 
if we have two boundaries together, there is indeed a SRE 
symmetric state on the two boundaries. The state is a 
product of entangled pairs of effective spins connected 
by a dashed line. The entangled pair can be chosen as 
1 00) + 1 11). In contrast to the single boundary case, we 
can locally project the two effective spins connected by 
a dashed line to the subspace |00) (00 1 + |11)(11| and on 
this subspace, the symmetry acts in an on-site fashion. 

This result should be expected because if we have two 
pieces of sheet with boundary and glue them back into 
a surface without boundary, we should have the origi- 
nal SRE 2D state back. Indeed if we map the effective 
spins back to the original degrees of freedom |0) —> 1 00) 
and |1) — > 1 11), we see that the SRE state between two 
boundaries is just the a chain of plaquettes |0000) + |1111) 
in the original state. 

This model serves as an example of non-trivial SPT or- 
der in 2D SRE states that only needs to be protected by 
on-site symmetry. In order to prove the special property 
on the boundary of CZX model and have a more com- 
plete understanding of possible SPT orders in 2D SRE 
states with on-site symmetry, we are going to introduce 
a mathematical tool called Matrix Product Unitary Op- 
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erator. We will show that 2D SPT phases are related to 
elements in H 3 (G, 17(1)) which emerge in the transforma- 
tion structure of the matrix product unitary operators. 
The definition of matrix product unitary operator and 
some basic properties are given in appendix D. The dis- 
cussion in the next section is general, but we will work 
out the CZX example explicitly for illustration. 

V. MATRIX PRODUCT UNITARY 
OPERATORS AND ITS RELATION TO 3 
COCYCLE 

In this section, we discuss the matrix product unitary 
operator (MPUO) formalism and show how the effective 
symmetry operation on the boundary of CZX model can 
be expressed as MPUO. Moreover, we are going to relate 
MPUO of a symmetry group to the 3-cocycle of the group 
and in particular, we are going to show that the CZX 
model corresponds to a nontrivial 3-cocycle of the Z 2 
group. 

A matrix product operator acting on a ID system is 
given by, 44 

0= Tr{T il ^T t2 ' l ' 2 ...T iN ^)\i' 1 i / 2 ...i' N )(i 1 i 2 ...i N \ 

(12) 

where for fixed i and i', T 1 ' 1 is a matrix with index a 
and j3. Here we want to use this formalism to study 
symmetry transformations, therefore we restrict O to be 
a unitary operator U. Using matrix product represen- 
tation, U does not have to be an on-site symmetry. U 
is represented by a rank-four tensor T£ q on each site, 
where i and i' are input and output physical indices and 
a, j3 are inner indices. Basic properties of matrix product 
unitary operators are given in appendix D. 

In particular, the symmetry operator Uczx (we omit 
the ~ label for effective spins in following discussions) on 
the boundary of the CZX model can be represented by 
tensors 

T^{CZX) = |0)(+|, 

T^{CZX) = |1)<-|, (13) 

other terms are zero 

where |+) = |0) + |1> and |-> = |0) - |1). It is 
easy to check that this tensor indeed gives Uczx — 
C Zi2...CZ niXi...Xn . 

The other element in the Z 2 group-thc identity 
operation-can also be represented as MPUO with ten- 
sors 

T°>°(J) = 10><0], 

T^(I) = \0){0\, (14) 
other terms are zero 

These two tensors are both in the canonical form as de- 
fined in appendix D. 

If two MPUO T(g 2 ) and T(g x ) are applied subse- 
quently, their combined action should be equivalent to 



T(gig 2 ). However, the tensor T(gi,g 2 ) obtained by con- 
tracting the output physical index of T(g 2 ) with the in- 
put physical index of T(gi), see Fig. 5, is usually more 
redundant than T{g\g 2 ) and might not be in the canon- 
ical form. It can only be reduced to T(gig 2 ) if certain 
projection P gi g2 is applied to the inner indices (see Fig. 
5). 





T(gi) 






p i 


r s v s 2 


T(gig 2 ) 




T(g 2 ) 



FIG. 5. Reduce combination of T((?2) and T(gi) into T(gi<?2). 

P gi) g 2 is only defined up to an arbitrary phase fac- 
tor e l6 ( 9l ' 32 \ If the projection operator on the right side 
Pgi,92 i s changed by the phase factor e iet ^ 9l,92 \ the projec- 
tion operator Pt g2 on the left side is changed by phase 

factor e~ ie(9l - 92 \ Therefore the total action of P gi . g2 
and Pg lt92 on T(gi,g 2 ) does not change and the reduc- 
tion procedure illustrated in Fig. 5 still works. Moreover, 
from the discussion in the appendix D, we know that this 
is the only degree of freedom in P git92 ■ Up to a phase fac- 
tor, P gi ,g 2 is unique (on the unique block in the canonical 
form of T{g 1 ,g 2 )). 

Let us illustrate how the reduction is done for the 
symmetry group (I, Uczx)- For example, if we apply 
UczxUczx the totally action should be equivalent to /. 
However the tensor T(CZX,CZX) is given by 

T°-°(CZX, CZX) = |01)(+- |, 
T 1 ' 1 {CZX,CZX) = \10)(- + \, (15) 
other terms are zero 

This tensor is reduced to T(I) if projection 

Pczx,czx = (|01) - |10))<0| (16) 

and its Hermitian conjugate are applied to the right and 
left of T{CZX, CZX) respectively. 45 Adding an arbi- 
trary phase factor e W(CZX,CZX) to p czx czx does not 
affect the reduction at all. By writing Pczx.czx in the 
above form, we have made a particular choice of phase. 

Below we list the (right) projection operators for all 
possible combinations of g\ and <?2 of this Z 2 group. 

Pi. i = |00)(0| 

Pczxj = |00)(0| + |10)(1| , , 

Pi,czx = |00)(0| + |10)(1| [U) 

Pczx,czx = (|01)-|10))(0| 

Note that in giving Pg 1 .g 2 we have picked a particular 
choice of phase factor e l0 ( 9l ' 92 \ In general, any phase 
factor is allowed. 

Nontrivial phase factors appear when we consider the 
combination of three MPUO's. See Fig. 6. 
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FIG. 6. Different ways to reduce combination of T(gz), T(g2) 
and T(gi) into T(<?i<?2i?3). Only the right projection operators 
are shown. Their combined actions differ by a phase factor 
0(51,52,33). 



There are two different ways to reduce the tensors. We 
can either first reduce the combination of T(g\), T{g 2 ) 
and then combine T{g%) or first reduce the combination 
of T{g 2 ),T{gj) and then combine T(gi). The two differ- 
ent ways should be equivalent. More specifically, they 
should be the same up to phase on the unique block of 
Tgi, 52,53- Denote the projection onto the unique block 



of T(g 1 ,g 2 ,g 3 ) as Q g 



We find that 



^91,32, 33(^3 ® Pgi,gi)P. 9132,93 = 
,32,93 



91,929.3 



(18) 



From this we see that the reduction procedure is asso- 
ciative up to a phase factor <f>(g\, 32, 53)- According to 
the definition of cocycles in appendix B, we see that 
0(51,52,33) forms a 3-cocycle of group G. That is, 
0(31,32,53) satisfies 



VI. NONTRIVIAL 3-COCYCLE OF MPUO AND 
NONEXISTENCE OF SRE SYMMETRIC STATE 

In this section we will show that a symmetry defined 
by a MPUO on a ID chain can have a SRE symmetric 
state only if the MPUO corresponds to a trivial 3-cocycle. 
Therefore, the boundary of the CZX model must be gap- 
less or have symmetry breaking. For this proof, we will 
be using the matrix product state representation of SRE 
states. 

Suppose that the symmetry on a ID chain is repre- 
sented by tensors T^ l g(g). WLOG, T(g) is single-blocked 
and in the canonical form as defined in appendix D. As- 
sume that it has a SRE symmetric state represented by 
matrices A\ which is also single-blocked and in the 
canonical form. For a review of matrix product state 
formalism including its canonical form and single-block 
property see appendix C. 

Based on the result in Ref. 46 and Ref. 47 we can show 
that (see appendix D) 



A 1 = V\^2t 1 ' 1 ' (g)A L ')V 



(21) 



where V'V = I and V is unique on the single block of 
Y^i' T 1 ' 1 (g)A l up to phase. This is saying that we can 
reduce the MPS obtained from £\, T 1 - 1 ' (g)A l> back to 
the original form A 1 by applying V* and V to the left 
and right of the matrices respectively. See Fig. 7. 





T(g) 








V 

A 




A 



0(32,33,34)0(31,3233,34)0(31,32,33) 



(19) 



Let's calculate 0(31 
erated by U C zx- 

0(7,7,7) = 1 
0(7, CZX, I) = 1 
0(7, CZX, CZX) : 
(j}\cZX, I, CZX) : 



0(3i32, 33, 34)0(3i, 32, 3334) 

32, 33) explicitly for the group gen- 



0(7, 7, CZX) = 1 
cp(CZX, 1,1) = 1 

4>{czx, czx, i) = 1 

cj){CZX, CZX, CZX) 



-1 

(20) 

We can check that is indeed a 3-cocycle. The last 
term shows a nontrivial — 1. This minus one cannot be 
removed by redefining the phase of P gi , g2 in any way. 
Therefore corresponds to a nontrivial 3-cocycle for the 
Z 2 group. 

What does this nontrivial mathematical structure im- 
ply about the physics of the CZX model? In the next 
section we are going to answer this question by prov- 
ing that MPUO related to a nontrivial 3-cocycle cannot 
have a short range entangled symmetric state. That is, 
the boundary of the CZX model cannot have a gapped 
symmetric ground state. It cither breaks the symmetry 
or is gapless. 



FIG. 7. Reduction of the combination of T(g) and A into 
A. Here T 1 ' 1 (g) is a MPUO, A 1 is a matrix product state 
symmetric under T l,z (g). 

For a fixed representation of the SRE state A 1 and fixed 
representation of the MPUO symmetry T(g), V is fixed 
up to phase. We can pick a particular choice of phase for 
V. 

Now we consider the combined operation of T(g\) and 
T(g 2 ) on A. See Fig.8. 
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B 1 ,g 2 j- 



} 



FIG. 8. Two ways of reducing the combination of T(<?2), T(gi) 
and A into A. Only the right projection operators are shown. 
Their combined actions differ by a phase factor v?(<?i,<?2). 

We can either first combine T(g 2 ) and A and then 
combine T(g\) and A or first combine T(g\) and T{g 2 ) 
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and then combine T(gig 2 ) and A. The right projection 
operator for these two methods differ by a phase factor 
<p(5i, 92)- This phase factor can be arbitrarily changed by 
changing the phase of P fll , 92 . For following discussions, 
we fix the phase of P gi , g2 and hence f(gi,g2)- 

This is all the freedom we can have. If we are to com- 
bine three or more T's with A, different reduction meth- 
ods differ by a phase factor but the phase factor are all 
determined by ip(gi,g2). Consider the situation in Fig. 
9, where we are to combine T{g^), T(g 2 ) and T{g{) with 
A. 



T(.g 2 ) Bl ' 
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T(g 3 ) 
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T(g> e Tl 






J v 123 

A 




A v 4 




(d) 


(e) 





T(gi) 




T(g 2 ) V 




T(g,) V 2 L 




A 4-1 


(c) 



■} 



3-cocycle of the Z 2 group, the system with boundary can- 
not have a gapped symmetric ground state. This shows 
that the CZX model has nontrivial SPT order protected 
by on-site Z 2 symmetry as we have promised in section 
III. 



VII. GENERALIZATION TO FERMION 
SYSTEM 

Due to the interest in fcrmion SPT orders in interact- 
ing systems in two and higher dimensions, in this section 
we are going to give a fermionic version of the CZX model 
which also has nontrivial SPT order protected only by 
on-site Z 2 symmetry. 

In constructing this model, first we identify each spin 
in the CZX model with a fermionic mode and the spin 
|0) state with the zero fermion state, the spin |1) state 
with the one fermion state. Each site then contains four 
modes (see Fig. 3). Denote the creation and annihilation 
operator on each mode as c] and Cj. 

Fermion system has an intrinsic fermion parity sym- 
metry which is an on-site Z 2 symmetry given by 



pf = na - 2ci« 



(25) 



FIG. 9. Different ways of reducing the combination of T(gs), 
T(g2), T(gi) and A into A. Only the right projection oper- 
ators are shown. Their combined actions differ by a phase 
factor written on the arrow. 



This Z2 symmetry is always preserved. 

Similar to the CZX model we define another on-site Z2 
f 

symmetry Uq ZX , which is going to protect the nontrivial 
SPT order. 



To change the reduction procedure in Fig. 9 (a) to that 
in Fig. 9(c), we can either go through step (b) or steps (d) 
and (e). If we go through step (b), the phase difference 
in the right projection operators is 



<P l {9i92,9z)'P 1 (9i,92) 



(22) 



On the other hand, if we go through steps (d) and (e), 
the phase difference in the right projection operators is 



<%i, 92, 93 )v 1 (9i,9293)f 1 (ff2,93) 



(23) 



But these two procedures should be equivalent as the 
initial and final configurations are the same whose phases 
have been fixed previously. Therefore, we find that 



0(51,52,53) 



y(gi, 5253)1(52, 53) 

f(9i92, 93)^(91, 92) 



(24) 



and 0(51, g 2 , 53) must be a trivial 3-cocycle (see Eq. 
B13). 

This finishes the proof that: A ID system with symme- 
try defined by matrix product unitary operators can have 
a gapped symmetric ground state only if the matrix prod- 
uct unitary operator corresponds to a trivial 3-cocycle. 

Because we have shown that the symmetry on the 
boundary of the CZX model corresponds to a nontrivial 



Tjf _ Tjf Tjf 

u czx — u x u cz 



where 



U 



f _ 
x " 



n( c i+ c *) 



(26) 



(27) 



is a particle-hole transformation and 



(28) 



1=1 



It can be checked that U x and U cz commute with each 



other and they both commute with Pf. Therefore U ( 



czx 



1 

commutes with Pf. U czx generates an on-site Z 2 sym- 
metry. 

The Hamiltonian of the system is again a sum of local 
terms around each plaqucttc. = ^2 H£. . 



Hi = -X{ ® F%' f ® P 2 d J ® P l 2 J ® Pl J 



(29) 



(see Fig.3(c))where X( acts on the four modes in the 
middle plaquette as 



x{ = C4C3C2C1 + 4444 



(30) 
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and P 2 acts on the two modes in every half plaquette as 

P 2 — c. l c\c l+ icl +1 + c\cic\ +1 c i+ i (31) 

For the remaining four modes at the corner, H?. acts as 
identity on them. It can be checked that the Hamiltonian 
satisfies the fermion parity symmetry and the on-site Z 2 
symmetry generated by U[, zx . Moreover terms around 
different plaqucttes commute with each other. 

The ground state is then a product of plaquette states 

l^ 4 > = (1+4444)10) (32) 

where |fi) is vacuum state on the four modes 1 ~ 4 
around a plaquette. The ground state is short range en- 
tangled and symmetric under both Pf and U^ zx . If 
Uq ZX can be violated, we can disentangle this state into 
a product of states on each mode without violating the 
fermion parity symmetry. However, if Uq ZX is preserved, 
this state is inequivalent from a trivial product state. 

The nontrivial-ness of this model can be seen again 
from the boundary. The boundary of this fermion model 
is the same as that of the spin model, because the ef- 
fective degrees of freedom in each half plaquette has two 
states: the vacuum state on the two modes and the fully 
occupied state on the two modes. These two states are 
both bosonic, therefore the boundary can be treated as a 
spin system just like for the CZX model. The symmetry 
on the boundary is again generated by Uczx which we 
have shown cannot have a SRE symmetric state. There- 
fore, this fermionic CZX model has a nontrivial on-site 
Z 2 symmetry protected SPT order. 

VIII. SUMMARY 

In this paper, we have given the explicit construction 
of a two dimensional interacting spin model with nontriv- 
ial on-site Z 2 symmetry protected topological order. We 
found that the system is highly nontrivial because if it 
has a boundary the boundary is either gapless or breaks 
symmetry. We showed this by writing the effective sym- 
metry transformation on the boundary as a matrix prod- 
uct unitary operator and revealed a nontrivial 3-cocycle 
structure in its transformation rule. We proved that any 
matrix product unitary operator related to a nontrivial 3- 
cocycle in H 3 (G, U(l)) cannot have a gapped short range 
entangled symmetric state. 

This model could have interesting implication for the 
study of topological phases using tensor network presen- 
tation. In the tensor network repretation of topological 
phases, it has been understood that in one dimension 
injective tensors provide a complete characterization of 
gapped ground states and its gauge transformation under 
symmetry reals the SPT order of the phase. 28 ' 29 In higher 
dimensions, similar analysis of injective tensors have been 
carried out. 48-50 However, the fact that the ground state 
wavefunction of CZX model has a loop structure and can- 
not be represented by an injective tensor tells us that we 



need to consider more general forms of tensors in order 
to study interesting SPT orders in more than one dimen- 
sion. Identifying the proper set of tensors for the char- 
acterization of gapped short range entangled phases in 
higher dimensions is an important open question. Or an 
alternative approach is to reduce the problem from 2D to 
ID by considering the tensor representation of effective 
symmetry transformation on the boundary, as was done 
in this paper. How the reduction can be done in more 
than two dimension is unknown. 

The ID boundary of the CZX model presents new chal- 
lenges to our understanding of ID systems. While it is 
a locally interacting system with Z 2 symmetry, it does 
not have a gapped symmetric phase like in transverse 
Ising model. Moreover, the gapless excitations cannot be 
gapped by breaking translational symmetry like in spin 
1 /2 chains. The peculiarity of this system originates from 
the fact that this ID system can only exist as the bound- 
ary of a 2D system and not on its own. Finding a proper 
field theory description of this system would expand our 
current understanding of ID physics. 

The relation between SPT order and cocycle is not 
accidental. Actually the pattern has shown up in lower 
dimensions. 29 In zero dimension, symmetric states are 
classified by ID representation of the group, that is, class 
of 1-cocycles in ^ 1 (G, (7(1)) and in one dimension SPT 
phases are classified by projective representations of the 
group, that is, class of 2-cocycles in H 2 (G,U(1)). Here 
we make a connection between 2D SPT order and class of 
3-cocycles in H 3 (G,U(1)). In fact, this relation is more 
general. In another paper, we are going to show that 
actually d dimensional SPT orders are related to (d + 
l)-cocycles in "H d+1 (G, C/(l)), which could lead to a full 
classification of SPT orders in any dimension. 

We would like to thank Zhenghan Wang for helpful 
discussions. This research is supported by NSF Grant 
No. DMR-1005541 and NSFC 11074140. 



Appendix A: Projective Representation 

Matrices u(g) form a projective representation of sym- 
metry group G if 

u{gi)u(g 2 ) = u(g 1 ,g 2 )u(g 1 g 2 ), 31,32 € G. (Al) 

Here uj(gi,g 2 ) e U(l) and uj(gi,g 2 ) 7^ 0, which is called 
the factor system of the projective representation. The 
factor system satisfies 

^{92,93)^(91,9293) = ^(91,92)1^(9192,93), (A2) 

for all 51,52,33 G G. If U)(g\,g 2 ) — 1, this reduces to the 
usual linear representation of G. 

A different choice of pre-factor for the representation 
matrices u'(g) = f3(g)u(g) will lead to a different factor 
system 0/(31,32): 

1 1 \ 0(9i92) , \ , Aq x 
w (31,32) = -r———uj{g 1) g 2 ). (A3) 
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We regard u'(g) and u(g) that differ only by a pre-factor 
as equivalent projective representations and the corre- 
sponding factor systems 0/(31,32) and u>(gi,g 2 ) as be- 
longing to the same class cj. 

Suppose that we have one projective representation 
u\(g) with factor system LJi(gi,g 2 ) of class lo\ and an- 
other 112(g) with factor system w 2 (gi, 32) of class u> 2 , ob- 
viously u\(g)®u 2 (g) is a projective presentation with fac- 
tor group uj\ (51, 32)^2(31, 52)- The corresponding class w 
can be written as a sum w\ + 0J2 • Under such an addi- 
tion rule, the equivalence classes of factor systems form 
an Abclian group, which is called the second cohomology 
group of G and denoted as H 2 [G,U(1)]. The identity 
element 1 € H 2 [G,U(1)} is the class that corresponds to 
the linear representation of the group. 

Appendix B: Group cohomology 

The above discussion on the factor system of a projec- 
tive representation can be generalized which give rise to 
a cohomology theory of group. In this section, we will 
briefly describe the group cohomology theory. 

For a group G, let M be a G-module, which is an 
abelian group (with multiplication operation) on which G 
acts compatibly with the multiplication operation (it the 
abelian group structure) on M: 

g-(ab) = (g-a)(g-b), g G G, a, b G M. (Bl) 

For the cases studied in this paper, M is simply the U (1) 
group and a an U(l) phase. The multiplication operation 
ab is the usual multiplication of the U(l) phases. The 
group action is trivial: 3 • a = a, g G G, a =G U(l). 

Let U! n (gx, g n ) be a function of n group elements 
whose value is in the G-module M. In other words, 
u> n : G™ -> M. Let C n (G,M) = {w„} be the space of 
all such functions. Note that C n (G,M) is an Abelian 
group under the function multiplication uj'^(gx, ...,g n ) = 
uj n (gi,...,g n )uj' n (gi,...,g n ). We define a map d n from 
C n [G,U(l)} to C n+1 [G,U(\)]: 

(d n u) n )(gi, ...,3„+i) = 

3i ' w„(3 2 , ...,3„ + i)cj r ( l _1) " +1 (3 1 ,...,3„)x 

71 

II^^ ( gi > —>9i-l>9i9i+l,9i+2>-9n+l) (B2) 
<=1 

Let 

B n (G, M) = {uj n \uj n = d n _ x a; n _i|a; n _i G C'-'fCM)} 

(B3) 

and 

Z"(G,M) = {u n \d n u n = l,uj n G C n (G,M)} (B4) 

B n (G, M) and Z n (G,M) are also Abelian groups which 
satisfy B n (G,M) c Z n (G,M) where B 1 (G,M) = {1}. 
The n-cocycle of G is defined as 

?T(G, M) = Z"(G, M)/B n (G, M) (B5) 



Let us discuss some examples. We choose M — U(l) 
and G acts trivially: g ■ a = a, 3 G G, a £ f (1)- In this 
case uj n (gi, ...,g n ) is just a phase factor. From 

(tii wi) (51,32) = ui (32)^1 (3i)/ w i (3152) (B6) 
we see that 

Z^G, C/(l)) = = ^(3132)}. (B7) 

In other words, Z 1 (G, U(l)) is the set formed by all the 
ID representations of G. Since B 1 (G,U(1)) = {1} is 
trival. H 1 (G,?7(1)) = Z 1 (G,U(1)) is also the set of all 
the ID representations of G. 
From 

(d 2 u 2 )(gi, 32,53) 

= ^2(52, 33)^2(31, 3233)M(ffi32, 33)^2(31, 52) (B8) 
we see that 

Z 2 (G,£/(1)) = { W2 | (B9) 
^2(32, 33)^2(31, 3233) = w 2 (3132, 33)^2 (31,32)}- 

and 

£ 2 (G, f7(l)) = {u 2 \u 2 (gi,g 2 ) = u>\ (32)^1 (3i)M (9m)}- 

(BIO) 

The 2-cocycle H 2 (G,U(1)) = Z 2 (G, J7(1))/£ 2 (G, C/(l)) 
classify the projective representations discussed in sec- 
tion A. 
From 

(^3X31,32,33,34) 

= ^3 (32 , 33 , 34 )^3 (31 , 3233 , 34 )^3 (31 , 32 , 33 ) , B11 s 
w 3 (31 32 , 33 , 34 ) W 3 (31 , 32 , 3334) 

we see that 

Z 3 (G,U(1)) = {uj 3 \ (B12) 

^3 (32 , 33 , 34 (3l , 3233 , 34 V3 (3l , 32 , 33 ) = X | 
w 3(3l32, 33, 34)^3(31, 32, 3334) ' ' 

and 

r3 /^ rm » r i / x ^2(32, 33)^2(31, 3233), 

6 J G,C/ 1 = {w 3 w 3 31,32,33 = —7 x — 7 r}, 

w 2 (3i32, 33)^2(31, 32) 

(B13) 

which give us the 3-cocycle W 3 (G,U(1)) = 
Z 3 (G,U(l))/B 3 (G,U(l)). 

Appendix C: Review: matrix product states and its 
canonical form 

In this section we review matrix product state and its 
canonical form which was first derived in Ref. 46. Similar 
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ideas are going to be used in the study of matrix product 
unitary operators. 

A matrix product representation of ID state is 



\i/>) = ]T Tr{A h A i2 ...A iN )\hi 



i2 — tN) 



(CI) 



Ai's are D x D matrices. 

Define double tensor E for the MPS as 



E = A i ® A * 



(C2) 



Equivalently, E can be expressed as a completely positive 
quantum channel £ as 



(C3) 



and the corresponding dual channel £* as 

£*(X)=J2A\XA i (C4) 

i 

The correspondence between E and £ , £* is as follows. 
Suppose that X and Y are D x D matrices which satisfy 



Y = S(X) 



(C5) 



Combine the two indices of the matrices into one and 
write them as vectors 

{Vx)(a-1)D+P = X a ,/3 (VV)(o-l)£>+^ = Y<*,P (G6) 

Vx and Vy are then related by E as 

w x = vv (C7) 

Similarly, if 

Y = £*(X) (C8) 

then 

(C9) 



V±E = V} 



We will use E and £, £* inter-changably, whichever is 
more convenient. 

From the structure of £ and £* we can put Aj's into 
a canonical form. Suppose that the largest magnitude 
of the eigenvalues of £ is Ai > 0. There could be mul- 
tiple eigenvalues \ie l6k of this magnitude. As shown in 
Ref. 51, e l9k form a group and they are the pth root of 
unity. To get rid of this, we can just group p sites to- 
gether and the eigenvalues of the largest magnitude will 
all be real and positive. We still label them as Ai. 

Because £ is a completely positive channel, at least one 
of corresponding fixed points A 



£(A) = AiA 



(CIO) 



is positive-semidefinite. Denote the support space of A 
as P. It can be shown that AiP = PAiP. 46 Decompose 



each Ai into four parts A t = PA^P + PAiP± + P±AiP + 
P^AiP^. P±AiP = 0. PAiP± may not be zero. How- 
ever, it does not contribute to the MPS, therefore we can 
remove it safely. After doing this, Ai is decomposed into 
two blocks and A is a full rank positive fixed point of 
£ P (X) = Y,i{PAiP)X{PAiPy with eigenvalue X v 
Because 

£ P (X) = Y,{PA i P)X{PA i P)^ = ^(A l P)X(A,P)t 

i i 

(Gil) 

every fixed point of £p (within space P) is also a fixed 
point of £ with the same eigenvalue. Therefore, Ai is 
also the largest eigenvalue of £p. Suppose that £p has 
another fixed point Z of eigenvalue Ai which is not pro- 
portional to A. WLOG, we can choose Z to be Hermi- 
tian. (This is because Y li ( A i p ) z (. A i p Y = there- 
fore J2i( A iP) Z H A iPY = A i^ f - And because Z is not 
proportional to A, at least one of the Hermitian matrices 
Z+Z^ or i(Z-Z') is not proportional to A.) Diagonalize 
the Hermitian matrix h~ 1 / 2 ZK~ 1 / 2 and get eigenvalues 
z\ > Z2 >■■■ It is easy to see that A — ^-2 is another non 
full rank positive fixed point of £p with eigenvalue Ai. 
Therefore we can repeat the previous process and turn 
PAiP into smaller blocks. 

Repeat this process for every block until (1) the chan- 
nel £p k of every block k has a largest positive eigenvalue 
A/j. There is a positive full rank fixed point Ap k within 
subspace (2) There is no other fixed point within Pk 
of the same eigenvalue. (3) The block P± = I — J2k Pk 
which does not have a positive fixed point for largest 
eigenvalue must only have zero eigenvalue. The block 
could be non-zero in general, but it does not contribute 
to MPS. Note that £ fc P k + P± = I, AiP k = P k AiP k . 
Written in the blocks P k and P±, Ai is upper(or lower) 
triangular. 

Now we look at each block k separately but from the 
dual channel perspective. We can similarly block di- 
agonalize A\ if non full rank positive fixed point ex- 
ists for the largest eigenvalue of £ Pk . For each sub- 
block projection P kjl , P k ,iA k t = P ktl A\P ktl . A\ can be 
turned into sub-blocks A k / = P kt iA^P k j. Note that, if 
Ap fc , = Pk,i^-p k P k ,i, 



Ei A 



k-.i 



A Phil ( A i' l Y = Af l A Pk (A^Y 

= P Kl AlK Pk {A^Pk,i 



(C12) 



Therefore, within each sub-block, Ap tl is still a positive 
full rank fixed point of £p k , with eigenvalue X k . As there 
cannot be positive fixed points of other eigenvalue, 
must be the largest. Similarly, if X k is a fixed point of 
£p k , P k ,iXkPk,i is a fixed point of £p kl with the same 
eigenvalue. 

Proceed similarly as for £, we can block diagonalize 
A k into A 1 *' 1 such that £p k i has only one fixed point for 
its largest eigenvalue which is full rank positive. 

Finally, we arrive at a canonical form, which is com- 
posed of blocks Pk and sub-blocks P k j. Within each 
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sub-block, the matrices satisfy (1) the channel £p k , has 
a largest positive eigenvalue. The corresponding fixed 
point is full rank positive. (2) There is no other fixed 
point within the sub-block of the same eigenvalue. (3) 
the dual channel £p k ; also has a largest positive eigen- 
value. The corresponding fixed point is full rank positive. 
(4) There is no other fixed point within the sub-block of 
the same eigenvalue. 

A generic matrix product state has only one block in its 
canonical form. 46 We will call these MPS single-blocked 
MPS. Single-blocked MPS represents gapped, short range 
correlated ID states. The single-block property is a 
generalization of the injectivity condition for MPS. 46 A 
single-blocked MPS is injective if the dimension of the 
matrices equals that in the canonical form. On the other 
hand, a single-blocked MPS might not be written in a 
canonical form. It is in general more redundant. To do 
the reduction, necessary steps involves projection onto 
the single block and re-labeling the basis. Any invertible 
operation within the projected space might be added. 
However, if the resulting canonical form is fixed, the re- 
duction operation is unique within the projected space 
up to an arbitrary phase factor. 

Appendix D: Matrix Product Unitary Operators 

Similarly to MPS, a matrix product representation of 
operators acting on a ID system is given by, 44 

0= Tr{T il ^T l2 ' l K..T iN ^)\i' 1 i / 2 ...i' N )(i 1 i 2 ...i N \ 

{ik},{i' k } 

(Dl) 

Here we restrict to unitary operators U as we want to 
discuss symmetry operations. Using matrix product rep- 
resentation, U does not have to be an on-site symmetry. 
U is represented by a rank-four tensor on each site, 
where i and i' are input and output physical indices and 
a, (3 are inner indices. 

Just like every matrix product state can be reduced to 
a canonical form. every matrix product operator can 
be reduced to a canonical form also. To do so, we just 
need to treat the two physical indices as one and apply 
the procedure described in appendix C. Similar to MPS, 
we can also define double tensor/ quantum channel for 
each matrix product operator. The double tensor of T is 

E = 5Z TM ' ® (T hi ')* (D2) 

The fact that T represents a unitary operator puts 
strong constraint on the form of T. WU = I ® ... (g> / is 
represented on each site by tensor 



could contain multiple blocks, but each block must rep- 
resent the same operator i" ® ... ® I and takes the form 
Afc<5i,i»|fe)(fc|. \k)(k\ is the projection onto the kth block, 
Xk is a number. Later we will impose further constraints 
on U to get rid of multi-block. 

First we want to show that we can write every MPUO 
in an single-blocked canonical form. That is, the canoni- 
cal form contains only one block. Suppose that we start 
with a canonical representation of the symmetry opera- 
tion. In general, the canonical representation could have 
multiple blocks. We are going to show that this is not 
necessary as different blocks represent the same unitary 
operation. 

Suppose that a canonical MPUO contains two blocks 

T" = Tft 8 T$ (D4) 

T m represents MPO Ox and Tpj represents MPO 2 (not 
necessarily unitary). U = 0\ + 2 . 

The corresponding T contains four blocks 

TP'*" = tU ' ® (T l ' v )* 

— T*'*' ffi T*'*' m T i,i ' eft T i,i ' (^5) 

- ^ [11] w 11 [12] W 11 [21] W U [22] 

Tnyfc/i represent MPO OkO\,. Because T represents I® 
I... <g) I, each of its block must also do. Therefore, 

OiO\ = O x O\ = 2 0\ = 2 0\ =1® I... ® I (D6) 

That is, 0\ and O2 represent the same unitary operator 
and there is no need for multiple blocks. In the following 
we will always assume that T is written in a canonical 
form with only one block. We will call this the single- 
block condition for MPUO. 

With the MPUO representation defined for each sym- 
metry operation, we now want to know how the repre- 
sentation changes when two or more operations are com- 
bined. 

First let's consider what happens when U is combined 
with U> , As we discussed before, this is represented by T 
which could contain multiple blocks \k5i,i" \k){k\ in the 
canonical form. Correspondingly, the double tensor of T 

E=j2 Ti ' i '®{ Ti ' i ')* = J2 Ti ' i ( D7 ) 

has multiple eigenvectors \k) with corresponding eigen- 
values Afe. 

Define the correlator between two sets of operator pairs 
{of ,of} and {oj.oj} to be 

(oi,o 2 )u = J2 mn Mo\ n oW5TdW r ) 

~(Em TrKWoTt/t))^ Tr(o2Uo%W)) 

(D8) 

On the one hand, written in terms of tensors, the corre- 
lator is expressed as 



T must be equivalent to S^f on each site. We can re- 
duce T to the canonical form. The canonical form of T 



{oi,o 2 )u = Tt(E..E[ 0i }..E[02}..E) 

- Ti{E..E[o 1 ]..E)Ti{E..E[o 2 ]..E) 
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where 

E M = E mt (°TY 2 ' l3 (5T) u ' n T n ' 12 <8 (T i4 ' i3 ) 
E lo 2 ] = E„i(°2) i2,i3 (5§) i4,il T il ' <2 (8 (T* 4 >* 3 )* 



(D10) 



This is the same form as the correlation function of oper- 
ators oi = E m o™ &> o™ and 02 = En °2 ® m a matrix 
product state with double tensor E. From our knowledge 
of MPS, we know that the correlator decays as ^/Ai)'. 

On the other hand, we consider for simplicity only 
unitaries U which preserve locality of operators exactly. 
That is, if o is supported on a finite number of sites, UoW 
is also supported on a finite number of sites, though the 
number may be larger. We do not consider the local op- 
erators with exponentially decaying tails. 52 Under this 
restriction, it follows that when {o™,6™} and {02,02} 
are far apart 



= E m „Tt«o?C7ofC7C/t 5 n C 7t) 

= Emn M(0?U5?W)) ® {0W0W ] )) 

= E mn Tr(d?U3rW)Tr{d5U5$Ui) 



(Dll) 



the correlator (01,02)11 must be zero if the separation 
is large enough. Therefore, A2 = 0. E has only one 
eigenvector and T has only one block in its canonical 
form. 

Now we want to use this property to show that the 
single-block condition is stable under combination of 
MPUO's. That is, if we start with two MPUO repre- 
sented by T a and T b with only one block in the canonical 
form, their combination 



rpC^il" V * rria,ii' rpb,i i 

1 aa\PP l ~ / j 1 <x.B L a'.fi< 



(D12) 



also has only one block in its canonical form. Of course, 
written as above, T c is not necessarily in the canonical 
form. Note that the discussion in the previous para- 
graphs is actually on the special case where T & ' M = 

In order to see this, we take the double tensor of T c 



E c = Ei,i" TC '"' ® (T c ' u ")* 

(D13) 

T a and T b both have one block in their canonical form. 
Denote the projection onto the blocks as P a and Pb- 



Pa = \r a& )(r^\ 
p b = 1^)^,1 



(D14) 



Being the only eigenvector of E a and E b , and 
\4>a'a') are positive full rank if written as matrices 
A b a , ~, . The only term that contributes to the trace of E c 
is 



This is also true for any power of E c . 

This special property of E c tells us that E c has only a 
single non-zero eigenvalue. Suppose E c = A o |0)(0| + ill, 
|0) is short for \ip a ) <8> \ip b ). Ty(E c ) = A . Moreover, 
Ti{E c ) k = A§. Because Tr(£; c ) fe = Ei( A *) fc , ^ can be 
shown that Ai = 0, Vi > 0. The fact that E c has a single 
eigenvalue in turn tells us that T° contains only one block 
in its canonical form, because otherwise, E c would have 
at least n 2 non-zero eigenvalues with n being the block 
number. 

Therefore, we have shown that if we start with the 
canonical single-blocked tensor representation of some 
unitary operators, the tensor obtained from their con- 
catenation still has only one block in its canonical form 
and is hence single-blocked. For single-blocked T we can 
always apply the procedure in Ref. 46 (also discussed in 
detail in appendix C) to reduce it to a canonical form. If 
we have multiple ways to do the reduction, they must be 
equivalent. More specifically, projected onto the unique 
block, the reduction operation is unique up to phase (if 
the final canonical form is fixed, not up to gauge). This 
phase factor is going to play an important role in our 
study of SPT orders. 

Similar reduction procedure applies when a matrix 
product unitary operator acts on a matrix product state. 
In particular, suppose T 1 ' 1 is a MPUO and A 1 represents 
a MPS which is symmetric under it. Suppose that T M 
and A 1 are both in the canonical form and have only one 
block. Because T M represents a symmetry of A % 



A 1 = J2 T 



A 



(D16) 



represent the same matrix product state as A . Moreover, 
because A is short range correlated and T % ' % does not 
increase correlation length, A is still short range corre- 
lated and it also contains one block in its canonical form. 
However, note that T" is a matrix and the inner dimen- 
sion of A is in general larger than that of A . Therefore 
A may no longer be in the canonical form. Some reduc- 
tion procedure needs to be done to bring A back to the 
canonical form. 

Suppose that P is the projection onto the single block 
in the canonical form of A 1 . Due to the uniqueness of 
the canonical form of a MPS, P must be of the same 
dimension as A 1 and PAP must be equivalent to A up 
an invertible transformation Q. ' That is 



A = QP^APQ- 



(D17) 



Denote V r = PQ- 1 and V x = QP, we get A = ViAV r . 
Moreover, ViV r — I, the identity on the inner dimensions 
of A . As Q is unique up to phase, Vi and V r are unique 
on the single block of A up to a conjugate phase factor. 
With slight abuse of notation, we will denote V r as V and 
Vi as and we have 



A 1 = V^AV 



{\r)®\^ h )m a \®^ b \) 



(D15) 



(D18) 
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